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The main challenge in computing inclusive cross sections and decay spectra in QCD is posed by kinematic thresholds. 
The threshold region is characterized by stringent phase-space constraints that are reflected in large perturbative 
corrections due to soft and coUinear radiation as well as large non-perturbative effects. Major progress in addressing 
this problem was made in recent years by Dressed Gluon Exponentiation (DGE), a formalism that combines Sudakov 
and renormalon resummation in moment space. DGE has proven effective in extending the range of applicability 
of perturbation theory well into the threshold region and in identifying the relevant non-perturbative corrections. 
Here we review the method from a general perspective using examples from deep inelastic structure functions, event- 
shape distributions, heavy-quark fragmentation and inclusive decay spectra. A special discussion is devoted to the 
applications of DGE to radiative and semileptonic B decays that have proven valuable for the interpretation of data 
from the B factories. 



1. Introduction 

The calculation of inclusive differential cross sections and decay spectra is amongst the most important and 
00 ■ well-developed applications of perturbative QCD. The main challenge in computing such distributions arises from 
\ kinematic thresholds, and specifically from the exclusive boundary of phase space. The threshold region is character- 
' ized by stringent constrains on real gluon emission, leading to large perturbative and non-perturbative corrections. 

This region is important for phenomenology in a wide range of applications. Well known examples are deep inelastic 
, structure functions in the limit where Bjorken x approaches 1 event- shap e distributions in e"'"e~ annihila- 

i-G I tion near the two-jet limit Drell-Yan 0,13, 0-0 or Higgs 00 

production in hadronic collisions 

^P^' near partonic threshold, or at small transverse momentum; heavy-quark fragmentation inclusive decay 

spectra |4lj|-|56j. In all these examples resummation and identification of the relevant non-perturbative corrections 
, can significantly increase the range of applicability of perturbation theory. 
^ ' Stringent constraints on real gluon emission are reflected in large Sudakov logarithms 57], namely perturbative 
corrections associated with mass singularities that become parametrically large near the exclusive limit. Owing to 
factorization, Sudakov logarithms can be resummed. It is generally the case, however, that upon approaching the 
\ threshold (while keeping the hard momentum scale Q fixed) non-perturbative effects become dominant. Sudakov 
resummation is therefore useful in a restricted range tliat is bounded from both ends: the logarithms need to be 
large enough to dominate, but in the close vicinity of the threshold, where the logarithms are indeed very large, the 
perturbative expansion as a whole breaks down, and it no longer approximates the physical distribution. 

At the perturbative level, a typical infrared and coUinear safe [a] differential cross section^ (or decay spectrum) 
takes the form 

1 da(Q, x) ^ ^(^^(q2))^(1 _ ^) ^ R{a,{Q^),x), (1) 



cTtot dx 

where x is some measured momentum fraction, the distribution has support for x < 1 and a; — 1 is the thresh 
old where the leading order distribution (5(1 — x) receives purely virtual corrections, H{as) = 1-1- vias + . . . 



^Although the problem addressed here is completely general, for simplicity we consider an infrared and coUinear safe distribution, 
where coUinear factorization is not needed. We further simplify the discussion assuming a single differential distribution where the hard 
kinematics is fixed, such as the thrust (T) distribution in e+e" annihilation l23ll where Q is the c ente r-of-mass energy and x = T (the 

threshold, x = 1, corresponds to the two— jet limit) or the photon— energy spectrum in B > Xs7 15.1 where Q = mt and x = 2Ej/mi, 

(the threshold, x = 1, corresponds to the maximal value of Ej for an on-shell b quark decay). 
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by 



ri{x)as 



contains terms of the form 



In (1 — a;)/(l — x) , where the plus distribution is defined 



dxf{x) 



ln'=(f 



x) 



1-x 



-I + 



dx{fix)-fil) 



hi'°(f-a:) 
1-x 



(2) 



where f{x) is a smooth test function. These distributions incorporate the cancellation between real and virtual mass 
singularities that are associated with the limit where a parton is fully resolved. At order a" Sudakov logarithms 
appear with < fc < 2n — 1. These perturbative terms get large as x — > 1 even if the coupling is small, and therefore 
resummation is necessary. 

The problem of the threshold region is best formulated in moment space, 



a(Q,7V) = 



dxx 



1 da{Q,x) 



(3) 



/o f^tot dx 

where the distributions become analytic functions of the moment index N . High spectral moments in Eq. Q are 
increasingly sensitive to the limit x — > 1 ; the plus distributions described above give rise to logarithmically-enhanced 



terms, In ' {N), which dominate at N 



Moment space is useful for Sudakov resummation because the multi- 



parton phase space factorizes there up to 0{\/N) corrections. Factorization properties of QCD matrix elements for 
soft and collinear radiation can then be used to prove exponentiation [57|-|62j: all the singular terms, a" In'^(A^), 
are generated by 



a{Q,N) = H{a,{Q'')) x Snd{Q,N) + 0{l/N), 



where the Sudakov factor takes the form 



^ oo n+1 V 

Sud(Q, TV) = exp i ^ ^ i?„,i a:{Q) In' {N) I , 



n=l 1=1 



(4) 



(5) 



where En^i are numerical coefficients, which depend on the relevant Sudakov anomalous dimensions and the coefficients 
of the 13 function. The resummed moment -space expression of Eq. can readily be matched to the fixed-order 
result to determine H{as{Q^)) and to account for the 0{1/N) contributions. Finally, the resummed cross section is 
obtained by an inverse Mellin transformation: 



1 da{Q,x) 



(Jtot 



dx 



2711 



a(Q,iV), 



(6) 



where the integration contour passes to the right of all the singularities of (t(Q. N). 

The physics underlying this resummation is the quantum-mechanical incoherence of dynamics at different momen- 
tum scales. This general principle is not restricted to the perturbative level — an observation that will be important 
for what follows. The typical excitations that are relevant at large N have virtualities of 0{Q^ /N) — the jet-mass 
scale or 0{Q/N) — the "soft scale". Because soft or collinear gluons do not resolve the hard interaction, they 
decouple from the process-dependent dynamics at the hard scale Q. 

For N — > oo the moments a{Q,N) factorize into "hard", "jet" and "soft" functions [1 B El El 0^- Upon 
introducing a factorization scale jip, the Sudakov factor in Eq. Q can be written as 



Sud{Q,N) = JN{Q;^iF) X SNiQ;nF), 



(7) 



where the jet function JNiQ'i (J-f) and the "soft" function Sn{Q; (j^f) sum up Sudakov logarithms associated with 
the jet-mass scale Q^/N and with the "soft scale" Q/N, respectively. The factorization-scale dependence^ cancels 



■^In the context of inclusive decays, see also Refs. |4.'l l4(l Is^l . 

^Factorization is implemented here using dimensional regularization (the MS scheme). In contrast with [It l34l . we do not employ strict 
scale separation; both J]v(Q;Mf) ^.nd SpfiQ', fJ-p) have some residual dependence on the hard scale Q; these functions are normalized 
such that their first moment is unity: Jjv=i(Q;A'f) = 1 and 5jv=i(Q;a*f) = 1- 
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In PT: on-shell b quark 

Figure 1: Factorization of inclusive decays into "hard", "jet" and "soft" functions, with virtualities of order mb, mb/^/N 
and nib/N, respectively. The jet function describes an unresolved jet with a given mass; see Ref. !^ and Sec. |21 below. In 
inclusive decays the "soft" function Sjv(jnt,; ^i?) has the process-independent interpretation of a quark distribution function. 
This function describes the longitudinal momentum distribution of an off-shell b quark field in the initial state (an on-shell b 
quark in perturbation theory, and a B meson in the full theory). See Ref. [s^ l for a precise definition. 



exactly in the product in Eq. Q . The phase-space origin of the logarithms is illustrated in Fig. ^ in the example of 
inclusive B decays. 

Sudakov logarithms, reflecting the hierarchy Q/N <C Q/Vn <^ Q, can be resummed to all orders. One expects, 
however, that a purely perturbative treatment will only be valid for Q/N ^ A. Even then, extremely soft gluons with 
virtualities C(A), whose interaction is not described by perturbation theory, have some effect on the distribution. 
These gluons become increasingly important as the hierarchy Q/N ^ A is removed. Considering the large-TV limit 
for fixed Q, one finds that the threshold region is characterized not only by Sudakov logarithms but also by non- 
perturbative effects that are inversely proportional to the "so/Z scale" , powers of NA/Q, which eventually always get 
larger than the logarithms as N increases. 

Importantly, the perturbative calculation itself reflects the presence of parametrically-enhanced power corrections: 
the series in the exponent of Eq. @ is non-summable owing to infrared renormalons In conventional 

Sudakov resummation this infrared sensitivity is ignored: the perturbative sum is truncated, not regularized. A 
regularization is necessary in order to systematically separate between logarithmically-enhanced perturbative cor- 
rections and power-enhanced non-perturbative power corrections. Such power corrections are uniquely associated 
with the phase-space region from which the logarithms arise, and therefore their physical interpretation is straight- 
forward. For example, in event-shape distribution the all-order resummation of (large-angle) soft gluon emission 
exposes parametrically-enhanced hadronization corrections; in inclusive decay spectra (Fig. ^ all-order soft-gluon 
resummation exposes parametrically-enhanced corrections distinguishing the quark distribution in the initial-state 
meson from that in an on-shell heavy quark. 

The presence of infrared renormalons and parametrically-enhanced power corrections in the moment-space Su- 
dakov factor is completely general, and it has far-reaching implications for the analysis of inclusive distributions. 
On general grounds, independently of the renormalon perspective, treatment of p ower corrections on the soft scale 
for A'^ ^ Q/A requires the introduction of a new nan- pe^t^a^e Jwact^on 0000, the so-called "shape 
function" . Unfortunately, the "shape function" is hard to constrain theoretically, while if constrained by data alone 
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the result would strongly depend on the functional form assumed and on the perturbative approximation with which 
it is combined. 

In practice, the threshold region was addressed by different tools in different applications, usually introducing some 
infrared cutoff on the perturbative result and parametrizing the contributions from the infrared as a "shape function" . 
Having little theoretical constraints on this function, the predictive power has been, in general, very limited. The 
most familiar examples are: 

• Event— shape distributions: Hadronization effects are taken into account though a leading non-perturbative 
"shift" of the Sudakov resummed spectrum 0,0,0 or, when addressing the peak region, through a con- 
volution with a "shape function" [isl l2lj . The general strategy of quantifying hadronization effects by power 
corrections represents an important advancement with respect to their description through fragmentation mod- 
els in event generators. This approach quickly led to successful phenomenology for average values of event 



svent generate 

Q-lii, 



shapes |2!|; but the challenge of describing the spectrum in the peak region could not easily be met. 

• Heavy— meson production cross sections: Fragmentation of the heavy quark into a meson is taken into 
account through a convolution with a leading-power fragmentation function |63L |6J| . Applied with or without 
Sudakov resummation, usually without any momentum cutoff, see e.g. j65| . 

• Inclusive B decay: The momentum distribution of the b quark in the B meson is taken into account through 
a convolution with a leading-power "shape function" |4l|, |43, l4,l[ 14(1 l47j . Applied with or without Sudakov 
resummation, using a momentum cutoff. 

Major progress in addressing this problem from first principals was made in recent years by Dressed Gluon Expo- 
nentiation (DGE) . In DGE the moment-space Sudakov exponent is computed as an all-order Borel sum, avoiding the 
usual logarithmic-accuracy truncation. Power-like separation between perturbative and non-perturbative corrections 
on the "soft scale" Q/N is achieved by the Principal Value prescription, avoiding any arbitrary momentum cutoff. 
This opens the way for making full use of the inherent infrared safety of the observable, which often extends beyond 
the logarithmic level. Infrared factorization and renormalization-group invariance strongly constrain the universal, 
all-order structure of the exponent. Additional, observable-dependent properties are derived from renormalon cal- 
culations of the corresponding Sudakov anomalous dimensions. These translate into constraints on the parametric 
form of power corrections that affect the exponentiation kernel. The end result is that the physical spectrum can be 
approximated by the regularized perturbative DGE spectrum, supplemented by a few power corrections of a known 
parametric form. This highly predictive framework has led to successful phenomenology in a variety of applications, 
includin g in particular event-shape distributions [2^ . heavy-quark fragmentation 39] and inclusive B decay 
spectra |53i l55l |. 

Our purpose here is to review the DGE approach from a general perspective. We will discuss the applications 
mainly to illustrate how the method works and distinguish between general and process-specific properties. We 
start, however, in Sec. 01 by shortly describing the motivation and the present status of the calculation of inclusive 
B decay spectra, where the application of DGE has proven useful for the interpretation of data from the B factories 
[so. 53. 54. 55,.j[56i]. In Sec. Olwe explain the foundation of the DGE approach using the example of an unresolved jet. 
This analysis |5j is directly relevant for the problem of deep inelastic structure functions at large Bjorken x, where 
the power corrections can also be viewed as resummation of the Operator Product Expansion (OPE) 0]. Next, in 
Sec. 01 we turn to discuss large-angle soft radiation. We explain how constraints on power corrections follow from 
the properties of the corresponding Sudakov anomalous dimension. The success of the resulting power-correction 
phenomenology is demonstrated using examples from event-shape distributions, heavy-quark fragmentation and B 
decay spectra. This is followed by a short summary of our conclusions in Sec. O 

2. Inclusive B decay spectra 

Amongst the most important contributions of the B factories are inclusive B decay measurements [gR IgtI : 
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photon ~ 



b quark 



s quark P 




jet 

m^=(P3 -qf 



Figure 2: The kinematics in B — > Xs'y: the photon, carrying momentum q (with — 0), recoils against the hadronic jet. 
In general, the virtuality of the b quark when decaying is a consequence of its "primordial" Fermi motion as well as of its 
subsequent interaction. In the on-shell approximation, where perturbation theory applies, the "primordial" Fermi motion, 
namely the interaction of the b quark with the light degrees of freedom in the meson is neglected; then the h quark momentum 
in the meson rest frame is p = {rrib, 0) while the momentum of the light degrees of freedom is fc = (A, 0), with A = Mb — rrib, 
where Mb is the meson mass and nib is the quark pole mass. 



• inclusive semileptonic decays, which provide the most accurate way to determine the CKM parameters |Vcb| 
and |V^,b|- The ratio IKib/l'cbl is one of the important constraints on the unitarity tria ngle |68l It is also 
an essential input in probing new physics, for example through B^-B'^ mixing, see e.g. 



• inclusive radiative decay, B — > Xg^, that directly probes the short-distance interaction responsible for flavor 
changing neutral currents. This decay, which occurs only through loops within the Standard Model, provides 



an important constraint on possible new physics 

scenarios HQ. 



Inclusive decays are theoretically favorable over exclusive ones owing to their low sensitivity to the hadronic structure 
of the initial and final states. 

Owing to irreducible backgrounds, experimental measurements of inclusive decays are limited to certain kinematic 
regions. This makes the theoretical calculation of inclusive decay spectra an essential ingredient on the way from the 
measurement to its interpretation. This is the case for both the radiative and the charmless semileptonic decays; the 
exception is B — > Xclv that can be measured over the whole kinematic range. In B — > Xg^, only hard photons 
with > 1.8 GeV, corresponding to Mx = \/{Pb - qY = \/Mb{Mb ~ 2E^) < 3 GeV, can be distinguished 
from the background; the measurement becomes statistically limited only for ^ 2.1 GeV. For the charmless 
semileptonic decay, B — > X^lv, the situation is yet more difficult: an upper cut on the hadronic invariant mass 
Mx < 1.7 GeV must be applied to avoid the overwhelming charm background. Thus, the determination of |Vub| 
from inclusive measurements strongly relies on the theoretical calculation of the spectrum. 

As dictated by the Born-level process, the typical hadronic momentum configuration in heavy-to-light decays is 
jet-like (see Fig.|2l: in the B rest frame the hadronic system has a high energy and a small mass, Mx ^ Mb- The 
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uT 



2.35 
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2.2 
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+ Belle data 
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Belle data 
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a,=0.222 




1.6 1.7 1.8 1.9 



Figure 3: B — > Xs^ theory vs. data: The average energy (left) and the variance (right) with a cut E~f > Eo, as calculated by 
DGE in [53I in the on-shell approximation (i.e. with no power corrections), varying the short distance parameters m"^ and 
a"^ within their error range. The result is compared with data from Belle 01 . Inner and total error bars show systematic 
and statistical plus systematic errors (added in quadrature), respectively. 



spectra peak near threshold^, making the understanding of this hmit absolutely essential. As discussed above, only 
the small-Mx region is experimentally accessible. 

An important observation underlying the theoretical description of inclusive heavy-to-light decay spectra, is that 
the decaying b quark is not on-shell |4lL l42l |43| . The Fermi motion of the b quark in the meson involves momenta 
of 0(A), and therefore one can expect 0{A) smearing of the perturbative spectrum by non-perturbative effects. The 
common lore, which developed based on this physical picture, is that computing the spectrum in the peak region is 
strictly beyond the limits of perturbative QCD. 

In fact, the situation is significantly better |5j0,EHE^- a systematic on-shell calculation involving resummation 
of the perturbative expansion yields a good approximation to the meson decay spectrum, one vifhich provides an 
excellent starting point for quantifying non-perturbative corrections. 

The on-shell approximation is physically natural because the heavy quark carries most of the momentum of the 
meson: the b quark virtuality, C'(A), is much smaller than the mass. When considering the total rate this translates 
into a systematic expansion in inverse powers of the mass, where the leading corrections are ©(A^/m^) and are 
numerically small. 

The next, crucial observation is that the on-shell decay spectrum is infrared and collinear safe, namely all its 
moments have finite coefficients to any order in perturbation theory. This means that non-perturbative effects, which 
make for the difference between the on-shell approximation and the physical meson decay, appear in moment space 
only through power corrections. Of course, when referring to the on-shell decay spectrum one must address the 
question of the summability of the expansion (as well as the precise definition of the on-shell mass), which brings 
about the issue of infrared renormalons. 

In general, the virtuality of the decaying b quark is a consequence of its "primordial" Fermi motion and of its 
subsequent interaction. This, however, does not preclude the resummed on-shell spectrum being a good approxima- 
tion to the meson decay spectrum. Obviously, the "primordial" Fermi motion and the "subsequent interaction" need 
to be defined. This definition amounts to a separation between non-perturbative and perturbative contributions, 
which can be done in a variety of ways. Cutoff-based separation procedures [zl| that are often applied, 

hinder the possibility of using of the inherent infrared safety of the on-shell decay spectrum. They necessarily asso- 



'For example, the B > Xs^ spectrum peaks near E.^ ~ mi,/2, see Fig. HI below. 
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Figure 4: The B — > Xs'y photon-energy spectrum (corresponding to the electromagnetic dipole operator O7) in the on- 
shell approximation as computed in fixed-order perturbation theory at NLO and at NNLO, and by resummed perturbation 
theory (DGE) 0, matched to NLO and NNLO, respectively. A common value of the (Principal Value) pole mass, 
rub = 4.88 ±0.05 GeV, has been used. This value was computed [53 based on the measured short-distance quark mass in MS, 
mp = 4.20 ± 0.04 GeV 0. 



elate a significant contribution with tlie non-perturbative "primordial" Fermi motion, which must be parametrized 
in a so-called "shape function" . Having little theoretical guidance, this parametrization is rather arbitrary, and the 
predictive power is very limited. The alternative used in DGE is separation at the level of powers, which is imple- 
mented using Principal Value Borel summation. In this way one can make full use of the inherent infrared safety of 
the on-shell decay spectrum, leaving only genuinely non-perturbative contributions, which distinguish between the 
(properly regularized) quark distribution in an on-shell quark and that in a meson, to be parametrized. Moreover, 
this parametrization is well guided by the theory: a small number of non-perturbative power corrections of a known 
form (see below) would suffice. In fact, with present experimental data no power corrections are needed — see Fig. 13 

It must be emphasized that the on-shell decay spectrum can only be an approximation to the physical spectrum 
when considered to all orders. In particular, both the on-shell mass, which sets the perturbative endpoint at any 
order in perturbation theory — > '7ib/2, or x = 2E^/m\, — > 1, and the spectral moments defined with respect to 
X as in Eq. ||2Jl, have a leading infrared renormalon at w = ^. The corresponding ambiguity cancels upon computing 
the sum of the series for rub and for the Sudakov factor using a systematic regularization — see Eq. (|30|l below. 
Examining the real-emission contribution to the spectrum in fixed-order perturbation theory (Fig. ^ one indeed 
observes huge corrections^ going from order to order^. Considering the resummed spectrum by DGE, matched to 
NLO and NNLO, respectively [fi^lF^ . one finds instead a remarkable stability, reflecting the fact that the dominant 
corrections have been resummed. Another striking difference between the resummed spectrum and the fixed-order 
one is their different support properties: the fixed-order result is a distribution that has support for < mb/2, 



^The 'no resummation' curves in Fig. |1] should not be considered too seriously for several reasons. First, the pole mass used at each 
order is the same, while in a strictly fixed-order treatment, it would b e differe nt (it diverges quickly). Second, the coupling is arbitrarily 
renormalized at mj,, while the typical gluon virtuality (the BLM scale iTdlTTD is significantly lower. 

®A full NNLO calculation of the normalized spectrum corresponding to the electromagnetic dipole operator was performed in Ref. l78l . 
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while the resummed resuh is a function that has support for < {mi, + 0{A))/2 ~ AIb/2, i.e. close to the physical 
spectrum. Evidently, resummation makes a qualitative difference. 

The ultimate test of the theoretical predictions is of course the comparison with experimental data. In B — > Xgj 
spectral data from Babar and Belle has become available shortly after the publication of the theocratical predictions 
by DGE 53]. Comparison of the first two spectral moments, computed with a varying lower cut > Eq, with 
Belle data is shown in Fig.|31 a similar comparison with BaBar data can be found in Ref. These results show 
that the resummed on-shell decay spectrum indeed provides a good approximation to the meson decay spectrum, 
as anticipated in Ref. Comparison of the computed B — > Xgj spectrum with data has several immediate 

applications: 

• Extrapolation of the partial width from the region of measurement to the whole of phase space, in order to 
confront it with the Standard Model prediction and provide constrains on physics beyond the Standard Model. 

• Precise determination of the b quark mass, see Ref. 67]. 

• Determination of the leading power corrections, the ones associated with the quark distribution in the meson. 
With the present accuracy, power corrections cannot be determined, but higher precision data are expected 
in the near future. Owing to their universality, these power corrections can readily be used in improving 
theoretical predictions for other inclusive decay spectra, notably the charmless semileptonic spectrum. 

A pressing issue, which is high on the agenda of the B factories is the precise determination of jKib]- As men- 
tioned above, inclusive measurements of charmless semileptonic decays provide the most promising avenue for this 
determination. The main obstacle is the large extrapolation from the region of measurement to the whole of phase 
space, which strongly relies on theoretical predictions for the spectrum. The strategy applied so far made direct 
use of the universality of the quark distribution in the meson, by first fitting an ansatz for the "shape function" to 
the B — > Xs"f decay data and then using it when computing the partial branching fraction of B — > X^lv within 
the region of measurement. As discussed above, in this approach there is little theoretical guidance on the functional 
form of the "shape function". Here the alternative presented by DGE is very attractive js^l Isfij] : in this framework 
the on-shell calculation, which depends only on the short-distance parameters, provides a good approximation to the 
spectrum. Moreover, prospects are high for improving this prediction further by higher-order (NNLO) calculations 
and by parametrization of the first few power corrections whose iV-dependence is known. 

The Heavy Flavor Averaging Group (HFAG) 67] has recently performed a first comprehensive study of jVub] 
from inclusive decay measurements using DGE. This analysis, which is summarized in Fig. [SI provides the most 
precise determination of JKib] so far. Since different measurements in Fig. El correspond to different kinematic cuts 
with extrapolation factors ranging between ~ 2 and 5, the consistency of the resulting values for jVub] provides 
an additional evidence that the underlying description of the spectrum is good. Direct comparison of the DGE 
calculation"^ of the triple differential B — > X^lv spectrum |^ (or its moments with varying cuts) with experimental 
data would be useful to constrain power corrections. Such analysis is complementary to the B — > Xs^ one in 
several ways owing to the potential Weak Annihilation contributions and the different (non-universal) subleading 
power corrections in the two processes [s^-li^, effects that are hard to quantify theoretically and that have so far 
been neglected. 

Having seen the advancement achieved by DGE in the calculation of inclusive decay spectra, let us now return to 
the theoretical foundation of the approach. Infrared and coUinear safe observables, such as inclusive decay spectra 
or event-shape distributions, typically involve two distinct sources of large corrections, the jet function and the soft 
fimction in Eq. 0. We begin in SecOlby considering the exarnple of the jet function Jn{Q, M_f), which is common 
to a large class of inclusive cross sections and decay spectra |22|, I^Sl In Sec. 01 we turn to the soft function 
<S'jv((9,A'-f)- 



'^The DGE calculation has been implemented numerically in C++ facilitating phase-space integration with a variety of cuts. The 
program is available at ,http : //www ■ hep ■ phy ■ cam ■ ac ■ iik/~andersen/BDK/B2U/ 1 
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Figure 5: |Vub| as extracted using DGE by the Heavy Flavor Averaging Group (HFAG) [gJ from inclusive semileptonic 
measurements by CLEO, Belle and BaBar. 



3. An unresolved jet and deep inelastic structure functions near the elastic limit 

Sudakov resummation is a manifestation of the infrared safety of inclusive distributions at the logarithmic level. 
However, the Sudakov evolution kernel conceals infrared sensitivity at the povi^er level, which only becomes ex- 
plicit once running-coupling effects are resummed to all orders. Consequently, Sudakov resummation involves power 
corrections that exponentiate along with the logarithms. 

To demonstrate this, consider first the jet function 0,J2l Jn{Q] I^f) that can be defined in a process independent 
way as the quark propagator in axial gauge (Eq. (2) in 5]). It obeys the following evolution equation: 



din Jq 1 ~ X 

J {os) is the corresponding scheme-invariant anomalous dimension. It is usually decomposed in the MS scheme 



as 



J ))=A )) + — = Cf^^ + ■■■, (9) 

where the conformal part, A (cks), is the universal cusp anomalous dimension, which is the coefficient of the singular 
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part, 1/(1 — in the splitting function 87]. A alone eontrols the factorization scale dependence of the jet function: 

^^^^^fi%^=^(a.(^|,))lniV. (10) 

The anomalous dimension B{as) ~ — |^ + . . . is associated with coUinear singularities El Q in the axial-gauge 
light-quark propagator, and its perturbative expansion is currently known to NNLO [s^ Isof . 

The jet function Jn{Q, (J-f) is also the Sudakov factor controlling the x — > 1 limit of deep inelastic structure 
functions 0, 

(Q2) = f dxx^'-^F^ix.Q^) ~ \H{Q;piF)JN{Q;f^F) + 0{1/N)] xq^ifip) + 0{NK^/Q^), (11) 



twist— two coefficient function 

where qNil-tp) is the quark distribution function, and H{Q; jip) is a hard function that does not depend on N . The 
evolution of the quark distribution function qnIpf) is controlled at large N by the cusp anomalous dimension A (as)'- 

'-^^^^ = -A{a^i,l))lnN + Oil) (12) 

such that the product in Eq. (|ll|l is invariant. 

The evolution equation of the structure function itself (at leading twist) is 



dlnFf (Q^) _ d\nJN{Q-^iF) d\nH{Q;^iF) 



Oil/N) (13) 



implying that log-enhanced terms in the leading-twist coefficient function of F2 exponentiate in moment space, to 
any order in perturbation theory. 

The cancellation between real (x^^^) and virtual (—1) contributions guaranties that the evolution kernel on the 
r.h.s of Eq. H13|l or Eq. ^ is finite at any order in perturbation theory, despite the 1/ (1 — x) singularity. Given 
that the anomalous dimension (as(/i2)) is known to a certain order (currently NNLO |88|. with a good estimate of 
the N^LO 89]) one can use the /3 function and integrate the r.h.s of Eq. © order by order, obtaining the following 
expansion: 



h.s{Eq. m} - cf EE^^^' (1^^)'' (14) 



j=i ;=i 

where the coefficients Cj^i with I — j at any order j (leading logarithms, LL) are determined by the LO term in J7 (ofg), 
those with I = j — 1 (next-to-leading logarithms, NLL) require also the NLO term in J'{as), etc. The standard 
approach to Sudakov resummation is based on expressing the sum in Eq. H14|) as an expansion with increasing 
logarithmic accuracy, namely 

r.h.s{Eq. ©} . ^ ±9.iX)[^y ; A. 1,^, (15) 

where f3o = ]^Ca — ^^f- Here^ 5'o('^) = ln(l ^ X), and higher-order terms gn{X) are determined based on the 
expansion of (ag) to N"LO. This sum is therefore truncated at n = rimax, the order at which {as) has been 
computed. 



*A well— known fact is that gn(A) have Landau singularities at A = 1, corresponding to Q'^/N = A^. Obviously, the perturbative 
analysis is valid only when the jet-mass scale Q'^/N is sufficiently large compared to A^ . H owever, when inverting the Mellin transform 
according to Eq. JHJ one needs to deal with this spurious singularity for any value of x 1901 . This problem is completely avoided in the 
DGE approach. 
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The next, crucial observation [5|,|2^|29j is that the sum in Eq. (|15|l does not converge. The same is true for any 
reorganization of the terms in Eq. H14|l because the coefficients Cj^i increase as j! at high orders. This is the effect 
of an infrared renormalon that dominates the evolution kernel in Eq. (jH)) at high orders. Importantly, , being 
an anomalous dimension, is expected to be free of any renormalon singularities^ . Thus, the one and only source 
of factorial divergence in Eq. ((SJ is the integration over the momentum fraction x near x — > 1. In this limit the 
coupling in (q^sI/^^)) is probed at extremely soft momentum scales — > 0. The non-existence of the sum in Eqs. 
(|15|1 or H14(l reflects sensitivity to the far infrared at the power level. 

A systematic way to quantify this infrared sensitivity is to regularize the divergence as a Borel sum. To this end, 
imagine that the anomalous dimension J [as) on the r.h.s. of Eq. ((SJ were known to all orders, so one could write 
its scheme-invariant Borel representation: 



(16) 



where T{u) is the Borel transform of the coupling^". Using Eq. H16|) in Eq. 
X integration arriving at: 



(or in Eq. Hl^fl) one can perform the 



rfln Jjv(Q;^f ) 
dlnQ2 



Cf 
/3o 



du 



A2 
Q 



^) Tiu)Bjiu) 



r(A^)r(-u) 
r(iv - u) 



(17) 



Indeed, one finds potential renormalon ambiguities at positive integer values of u arising from the x — > 1 limit. 
As mentioned above, Bj{u) is not expected to have any renormalon singularities, however, unless it vanishes at 
u = k where k = 1,2,..., the evolution kernel (the r.h.s. of Eq. H17|l ) will have power-like ambiguities that scale at 
large A'' as [Nh? /Q'^)^ . Upon solving the evolution equation for Jn{Q; f^p), any such power terms would obviously 
exponentiate along with the logarithms. 

The discussion of renormalons can be made concrete by focusing on the gauge-invariant set of radiative corrections 
corresponding to the large-/3o limit^^. In this limit one can get analytic results for the Borel transform. The result 
for the anomalous dimension reads j^: 



BAu) 



5, 

ea 



largo (3o 



Sm TTU 
TTU 



1 



1 



1 - u/2 



(18) 



implying that renormalon singularities in Eq. H17|l at u = 1 and u — 2 are indeed present"'^^. Consequently, Su- 
dakov resummation cannot be considered a purely perturbative issue: the evolution kernel itself is only defined in 
perturbation theory up to powers. Since the l.h.s of Eq. I|13|) is an observable, these ambiguities must cancel once 
non-perturbative effects are systematically included. Indeed, within the OPE renormalon ambiguities reflect the mix- 
ing under renormalization between operators of different twist 0, ; in the case of F2 this issue is understood in 
full detail at twist four, see Ref. Of course, the OPE cannot be directly applied to the evolution equation (|17|l . 
which was only derived at the leading-twist level. The exponentiation of these power corrections i5| goes beyond 
the OPE analysis: it amounts to resumming the OPE. 

Importantly, Bj(u) in Eq. H18|l vanishes at certain integer values u = k, for fc > 3, so the corresponding renormalon 
ambiguities in the evolution kernel H17I) are absent. Since a single dressed gluon is a natural approximation to the 
evolution kernel (e.g. iteration of a single chain is generated by solving the equation, through exponentiation) we 
expect that the pattern of renormalon singularities exposed in the large-/3o limit would not be modified. We will 



®This is a conjecture. Explicit calculation in the large-/3o limit (Eq. I18i below) supports it. 

'^"In general we use the scheme-invariant Borel transform i91il where T{u) is the Borel transform of the two-loop coupling, see Eq. 
(2.18) in Ref. In the large-/3o limit T{u) = 1. 

^^Results in the large-/3o limit are obtained |92| by first considering the large-A^j^ limit, in which a gluon is dressed by any number of 
fermion-loop insertions, and then making the formal substitution Nf > — 6/3o. ^^^^ 

^^Eq. 1181 is consistent with renormalon calculations of deep inelastic coefficient functions, done independently of the x ^ 1 limit I94ll95l . 
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therefore assume that the zeros of Bj{u) in Eq. IjlSI) are the zeros of this function in the fuU theory. As we shall see 
in the next section, other Sudakov anomalous dimensions are characterized by their own pattern of zeros. 

In general, Borel singularities in QCD differ from their large-/3o limit in both the value of the residue and the nature 
of the singularity 92, 96] . However, Borel singularities in the Sudakov evolution kernel rcmain^"^ simple poles^'^ 
This follows directly from Eq. I|17|l and the above conjecture that Bj{u) itself does not have Borel singularities. 

To make practical use of Eq. H17() one obviously needs sufficient knowledge of the Borel transform of the Sudakov 
anomalous dimension, Bjiu). Perturbation theory directly gives the Taylor expansion of Bj{u) near the origin. 
Examining the integral in Eq. H17|l one finds that at small N only the immediate vicinity of the origin is relevant, 
while for N ^ 0{Q^ / A^) contributions from u ^ 1 and u ^ 2 become important. Owing to the strong suppression 
of contributions from large u through r(— u), the behavior of Bj{u) at large u is irrelevant, even for high moments. 
Aiming at power accuracy, one would need to have Bj(u) under control at least for < u ^ 1, and specifically, 
constrain the value oi Bj(u — 1) that controls the magnitude of the leading renormalon residue in the evolution 
kernel. 

It is natural to express Bj{u) in terms of its known large-/3o limit as foUows^^: 

Bj{u) = i?j(«)li„gc/3o X ; = 1 + '^("//Jo), (19) 

where Vj{u) accounts for contributions that are subleading for [3q — > oo; its Taylor expansion around it = can be 
determined order by order in perturbation theory. Eq. (|19() is consistent with B j {u) having the same pattern of zeros 
as exposed by the large-/3o limit 1)18(1 provided that Vj{u) itself has no singularities nor zeros, at least for positive 
integer values of u. A concrete ansatz^^ that satisfies these requirements is Vj{u) — exp(— wiu — W2U^ + •••)• Here 
the coefficients Wi are determined order by order, up to the order at which J^{as) has been computed 1^: based on 
the NLO result for the cusp anomalous dimension and Eq. (|5J) one finds wi = ^ (^j^ — , from the NNLO result 
for J{as) one can determine W2, see Eq. (2.33) in [s^, and finally from the approximate N^LO coefficient for 
J^{as) one can determine to good accuracy. 

The characteristic feature of the power ambiguities of Eq. (|17|l is their enhancement by the moment index iV, 
reflecting the fact that the physical scale involved is Q^/iV, the invariant mass of the jet. In the language of non-local 
lightcone operators [9^ , the large- A'^ limit corresponds to the limit of large lightcone separation between the quark 
fields 4]. Alternatively, within the local OPE this parametric enhancement is understood as power increase in the 
number of relevant local operators. This leads to an explosion in the number of non-perturbative parameters, making 
the description of structure functions in the large-x limit within the framework of the OPE completely impractical. 

Physically, the breakdown of the OPE and the appe arance of non-perturbative corrections on the jet mass scale 
are related to the transition into the resonance region lOol lOlj . where eventually the jet becomes a few exclusive 



^^Note that these are simple poles only in the scheme-invariant formulation of the Borel transform, where the Borel variable is conjugate 
to ln((3^/A^). In the standard formulation, where the Borel variable is conjugate to 7r/(/3o«s) the singularities turn into cuts, controlled 
by /3i//3g. See Eq. (B.13) in Ref. Isl . 

^*It implies I3|, in particular, that the mixing of higher-twist operators into the leading twist involves at large N pure powers, with 
no additional logarithms. In other words, one expects that the effective anomalous dimensions of higher— twist operators would have the 
same large- A'^ asymptotic behavior as the leading— twist operator, Eq. 1121 . 

Anomalous dimensions of higher-twist operators can of course be computed directly, independently of their power-like mixing with the 
leading twist. A significant progress in addressing this problem was made in recent years using the integrability property of lightcone 
operators l9Sll . A general result is that the spectrum of higher-twist anomalous dimension s, co mputed as a function of N , is bounded from 
below by the leading-twist anomalous dimension — see e.g. Eqs. (4.4) and (5.3) in Ref. l99l . This implies that the above expectation is 
indeed realized asymptotically. 

^^Starting with the NNLO result for Bj(u), Ref. Q considered several different functional forms that are consistent with Eq. 1181 . 
concluding that a multiplicative form as in Eq. I19i . where Vj{u) does not introduce any Borel singularities, results in good apparent 
convergence of the logarithm ic a ccuracy expansion, whereas Eq. 1181 with additive O{l/l3o) corrections would lead to extremely large 
subleading logarithms. Ref. l89l has show n t hat the logarithmic accuracy expansion converges well. This supports the multiplicative 
form of Eq. I19i . The new N^LO results l89l constrain Bj{u) further, and they may well facilitate an approximate determination of 
Bjiu = 1) corresponding to the leading renormalon residue in Eq. I17i . 

^^When using Eq. 1191 for phenomenology, the sensitivity of the result to the particular ansatz for Vj{u) represents the theoretical 
uncertainty due to unknown higher-order corrections. 
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states. Obviously, the wealth of non-perturbative phenomenon in this region cannot be described in detail using the 
OPE. However, parton-hadron duality may well work in moment space up to very large N . 

The discussion above suggests a clear track along which the evolution of structure functions at large-x should be 
computed 

• First, the evolution kernel on the r.h.s of Eq. IjlTI) needs to be regularized, e.g. by taking the Principal Value 
of the Borel sum. Upon solving the evolution equation with the regularized kernel (and matching with the 
fixed-order result) one obtains an improved approximation to the evolution of the structure function at large iV, 
which is valid up to ©(iVA^/Q^) corrections. 

• For sufficiently large iV, the 0{N / Q"^) power ambiguities become relevant, and the perturbative evolution 
should be supplemented by non-perturbative power corrections. Since power ambiguities appear in the evolution 
kernel, power corrections appear in the Sudakov exponent. The perturbative jet function is then modified by a 
multiplicative non-perturbative factor of a particular form that is dictated by the ambiguities^^, 

where J^^^^^ (Q) = exp (-uo^ [N - 1) (AVQ^) _ ^pV(^ _ ^^(^ „ 3) (AVQ')') ■ (20) 

The non-perturbative parameters 0JY2 correspond to the regularization prescription (Principal Value) applied 
to the r.h.s of Eq. ((T7jl . such that the product J^{Q, ^_f) x J^^^^^^ [Q) is prescription-independent. It should 
be emphasized that by introducing power terms according to the ambiguity structure we explicitly assume a 
minimal model for the non-perturbative contribution: the first two power terms must definitely be there, but 
it is impossible to exclude by present theoretical tools higher power terms. To the extent that this assumption 
holds, Eg. JI20() sums up the dominant corrections in the OPE, scaling as powers of iVA^/Q^, to all orders (c/. 

Eq. dni) B"^ 

(Q2) ^ [H{Q■^iF)Jl^{Q■^JiF) X J^^P^) (Q) + 0{l/N)\ x g^(^f.) + 0(K^/Q^). (21) 

Ultimately, la rge- x structure function d ata can be an alyzed along these lines. Conventional structure function 
phenomenology |l02j (see, however, Refs. does not focus on constraining high moments, and parton- 

distribution-function fits are typically done with a cut on the hadronic mass, Q'^{\ — x)/x > 12.5 GeV^, restricting the 
data considered. In this way one avoids the need for Sudakov resummation or higher-twist corrections. Obviously, 
such restrictions limit the possibility of constraining the parton distribution functions. With the theoretical tools 
described above, there are high prospects for improving these constraints. 



4. Soft gluons, Sudakov anomalous dimensions and parametrically-enhanced power 
corrections 

In the previous section we discussed the DGE approach in the context of the exclusive limit of deep inelastic 
structure functions. This problem is special since it is characterized by a single hierarchy: the invariant mass 
of the jet <5^(1 — x)/x gets significantly smaller than the hard momentum transfer Q^; when the jet-mass scale 
approaches the hadronic scale, the structure function directly probes the exclusive limit, where the jet is composed 
of a few hadrons. As explained in the introduction, other inclusive distributions have a more complex structure. 



^''Note that although the analysis is formally valid up to 1/A'' terms, we keep the full A'' dependence of the ambiguities. This is a 
convenient choice; for example it implies that the normalization convention of {Q, fip) (vanishing at A'^ = 1) are then satisfied by 
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where the threshold region is characterized by a double hierarchy. Specifically, infrared and coUinear safe observables 
such as event-shape distributions and inclusive decay spectra involve the following three scales: hard 0{Q), "jet" 
0{Q / y/N) and "soft" 0{Q/N). In this situation the threshold region, where Sudakov resummation is important, is 
parametrically wider than the exclusive region discussed above. Moreover, the most important power corrections are 
then the ones associated with large-angle soft emission, powers of NA/Q, whereas the non-perturbative structure of 
the jet, i.e. powers of NA'^/Q^, can usually be neglected. 

In this section we consider Sudakov factors associated with soft-gluon emission, which we denote in general by 
Sn{Q', ^J'F)■ The Sudakov factor SNiQ',^J■F) is defined at the level of the cross section so it involves the amplitude 
times the complex conjugate amplitude. Soft gluons do not resolve the hard interaction nor the structure of the 
jet; they couple to the color field along the classical trajectories of hard partons P, IbcI 0|. The relevant log- 
enhanced term s can therefore be co mputed in the Eikonal approximation or, equivalently, by considering Wilson-line 
operators [s^, lOSl 107 . llOSl Il09| — see Fig. Such field-theoretic definiti ons are useful because they readily 



generalize the perturbative Sudakov factor to the non-pertubative level |18[ l52l llf 

"Soft" Sudakov factors corresponding to different observables are the same only in the conformal limit (or to 
leading logarithmic accuracy) where just the cusps along the Wilson loop contribute; subleaing logarithms and power 
corrections depend on the geometry of the process (Fig. and on the way particle momenta are weighted. Indeed, 
already at the perturbative level Sn{Q',IJ'f) have different physical content in different processes, for example: 

• In Drell-Yan or Higgs production near partonic^^ threshold Sn{Q', IJ^f) describes the energy distribution of soft 
gluons (1 — a; is the small energy fraction carried by gluons) emitted by the incoming partons, which can be 
represented in terms of lightlike Wilson hues HSSElliSEiEiEi- 



• In a large class^^ of event-shape distributions in e+e" annihilation |22L l23l |2J|, S^iQ^fJ-F) stands for the 
contribution to the observable by (large-angle) soft gluons that are emitted off the recoiling quark-antiquark in 
the two-jet limit. It can formally be described in terms of a weighted integral over a matrix element involving 
lightlike Wilson hues [HEi- 

• In heavy-quark fragmentation SM{m] p,p) describes soft gluon radiation off an off-shell heavy quark that was 
initially produced with high energy of 0(Q), where Q ^ m, and then gradually looses its energy approaching 
its mass shell. This radiation is dominated by transverse momenta (or virtualities) of order \k±\ ~ rn(l — x), 
— the "soft" scale — where x is the longitudinal momentum fraction carried by the final -state on- s hell heavy 
quark (which represents the detected heavy hadron on the perturbative level) S Ei EiB, EH Eli The 
description of SN{m; jip) in terms of Wilso n lin es involves timelike lines connected by a finite lightlike segment 
of length that is proportional to N [s^ . llOTf . 

• In heavy-meson decay spectra S^i'm] hf) describes the soft gluon interaction with the heavy quark prior to 
its decay. Formally, SN{m] fip) is the Sudakov factor associated with the longitudinal momentum distribution 
of an off-shell heavy-quark field in an on-shell heavy-quark initial state. Remarkably, it is identical to the 
Sudakov factor in heavy-quark fragmentation , while the non-perturbative dynamics in the two cases is 
different. 

Despite the different physics described by the "soft" Sudakov factor in different processes, SN{Q',fJ-F) has a universal 
structure that is constrained by factorization, infrared safety and renormalization-group invariance. In particular, 
in all the examples mentioned above, upon neglecting 0{\/N) corrections one can establish an evolution equation 



'^^If the initial— state parton distribution functions are defined in dimentional— regularization based factorization schemes (such as MS), 
these cross sections involves only a "soft" Sudakov factor a jet-mass scale never appears. 

'^^We consider observables where contributions from individual soft gluon emissions are additive, leading to exponentiation in moment 
space. 
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event-shape 
weight 




Heavy-quark distribution (b decay) Heavy-quark fragmentaion (b production) 



Figure 6: Representation of the soft Sudakov factor Sn{Q', f^p) in terms of Wilson-line operators in Minkowski space-time, 
where the time is the horizontal axis. The four diagrams correspond to Drell-Yan or Higgs production, event-shape dis- 
tributions in the two-jet limit, heavy-quark distribution (inclusive decays), and heavy-quark fragmentation (heavy-meson 
production). In each case the lines close at infinity to form a loop, which guarantees gauge invariance. 



that is analogous to the jet-function case ©: 

This equation implies exponentiation of all the logarithms to any order in perturbation theory. The observable- 
dependent, but scheme-invariant anomalous dimension S (as) is usually decomposed in MS as: 

6 [as) = A[as) + - ^ — = Cf h . • . , (23) 

m /i^ dag TT 

where A{as) is the universal cusp anomalous dimension appearing in Eq. Q and V^ag) is observable-dependent. 
The dependence of Sn^Q] (J-f) on /x_f is the same as that of the MS quark distribution, i.e. it has the opposite sign 
to Eq. 1)10(1 above, making the product in Eq. (O factorization-scale invariant. 

Similarly to the jet-function case, infrared sensitivity at power level is concealed in the r.h.s to the evolution 
equation (|22|) : in the x — > 1 limit the coupling is probed at extremely soft momentum scales, fi^ — Q^{l~x)^ — > 0. 
In conventional Sudakov resummation this sensitivity is not regularized and the resulting divergent series is truncated 
at some order (some logarithmic accuracy). In order to regularize the divergence and quantify the infrared sensitivity 
as power-suppressed ambiguities we introduce a Borel representation^*^ for the anomalous dimension S [as) in analogy 



^"As before we extract the color factor out of the Borel transform. Cp in Eq. 1241 refers to radiation off hard quarks; For Higgs 
production though gluon-gluon fusion it should be replaced by Ca- This issue aside, Sn{Q; fJ-p) is identical to the Drell-Yan case. 
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with Eq. ((TB)l : 



5(a,(/i^))=^y^ dui^-j T{u)Bs{u). (24) 
Using H24|l in (|22|l . and performing the x integration, one obtains the DGE form of the evolution equation: 



dlnQ2 



f3o 



= / du[^] T{u) X Bs{u) X 



r(jv)r(^2^) 1 

r(iV-2u) 2u 



(25) 



Note that while Eq. p7|) indicates potential renormalon singularities at all mfeger values of w, in Eq. (|25|1 renormalons 
may occur also at half integer values. The corresponding power ambiguities at u = ^ {k being a positive integer) would 

i , which would result in 



scale at large N as {NK/Q)^ . Specifically, the leading renormalon ambiguity can be at u 
a substantial 0{NA/Q) effect — a global shift of the distribution. Conse quen tly, these power corrections received 
much attention in various applications including Drell-Yan production (291 . Isil l32l | (where the u = ^ ambi guit y was 
found not to occur, see below), heavy-quark fragmentation [sJIs^EI and especially event-shape variables j^-Q. 
Notably, the relevance of renormalons in computing inclusive B decay spectra has not received proper attention^^ 
until recently 0, El El El El El 

Since S (as) is an anomalous dimension, Bs(u) is not expected to have any renormalon singularities of its own. 
It follows that the sole origin of renormalons in the evolution kernel is the integration over x for x — > 1, leading 
to simple poles in Eq. (|25p ^. As usual, perturbative QCD calculations facilitate the determination of the first few 
orders in S (as), corresponding to the expansion of Bs{u) around the origin. In order to know the pattern of infrared 
renormalon singularities of the evolution kernel, one needs to know Bs(u) as an analytic function, and specifically 
find its zeros in order to deduce which power ambiguities are missing. Here the large-/?o limit becomes useful: when 
considering the evolution kernel a single dressed gluon becomes a natural approximation; non-Abelian diagrams are 
expected to modify the residue of the renormalons by Cyi//3o-dependent contributions (0(1) effects) but not to alter 
their positions. Thus, the pattern of zeros of Bs{u) can be directly deduced from the large-/3o limit. Table |l] lists 
the results in some important examples. For each observable we quote^^ the expression for bg{u) that is related to 
the corresponding Bs{u) by: 



Bs{u) 



e3' 



sm TTU 



bs{u) X Vs{u); 



Vsiu) = l + 0(u//3o), 



(26) 



where Vs{u) has no renormalon singularities nor zeros. The table also shows where, on the positive real axis, 
-^<s(''^)liargc A, vanishes, thus eliminating the potential renormalon singularity on the r.h.s. of Eq. (|25|) . Finally it shows 
which power ambiguities do appear. Evidently, the renormalon pattern is different for each class of distributions. 

In order to make practical use of the DGE formulation for the evolution equation H25|l . one needs to incorporate 
the known coefficients of the Sudakov anomalous dimension into an analytic function of u, namely write an ansatz 
for Ks(u) in Eq. if^ . A concrete example is Vs{u) = exp(— j/iu — y2u'^ + ■■■), where Vi = ^ — 5) is universal 
while Hi for i>2 are observable-dependent. These coefficient are determined order by order in perturbation theory 
up to the order at which S{as) has been computed: from the NNLO result for S{as) one can determine 2/2, etc. 
The sensitivity of the resulting Sudakov factor Sn{Q', f^p) to the particular ansatz for Vs{u) reflects the remaining 
theoretical uncertainty due to unknown higher-order corrections, see e.g. Eq. (3.27) in Ref. js^ and the discussion 
that follows. 



^'^See, however, earlier work in Ref. Ill3l . 

Going beyond the large-/3o limit the singularities remain simple poles only upon using the scheme— invariant formulation of the Borel 
transform 1251 . where the effect of the two- loop running coupling is factored out into T{u). In a standard Borel representation with 
respect to the coupling, th e s imple poles will be replaced by cuts whose strength is determined by the first two coefficients of the /3 
function, see Eq. (B.13) in l53l . 

■^^Technically, bs{u) is dir ectly obt ained from the Laplace transform of the "characteristic function", which is the leading— order r esul t 
with a single off-shell gluon Ill5llll3l . An example of such a calculation (for the C parameter in e+e~ annihilation) can be found in |24| . 
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Table I: Results for Sudakov anomalous dimensions in the large-/3o limit, Eq. ^ B El El Q S El El El, and the 
corresponding power corrections in the evolution kernel, the r.h.s of Eq. 1251 . 



Observable 



bs{u) Bsiu) = for It > power corrections 



Drell-Yan or Higgs production 
near partonic threshold 



Event Shapes (e+e — > jets) 
near the two-jet limit 



Heavy Jet Mass 
Thrust 
c parameter 



. r^i-«) 
'r(i -2m) 

1 

2 

. r^(i + ^) 
'r(i + 2u) 



1 3 
2' 2' 



M= 1,2,. 



— 1 , even 



^ I , odd 



Heavy Quark Fragmentation for x ~* 1 

Heavy Quark Distribution for x 1 
(the hard scale Q is the quark mass m) 



(1 



u = 1 



I I . K 7^ 2, mteger 



\ m J 



Examining the Borcl integral in Eq. (|25|l one finds that owing to the suppression by r(— 2u), the large u region 
is irrelevant. For N ^ Q/A the integral is dominated by small values of u, while for large N, of 0{Q/A), the 
intermediate region, u ~ i, and above becomes relevant. In this respect additional constraints on Bs{u) away from 
the origin are very useful in extending the perturbative treatment to larger N. Such constraints are provided for 
example by the zeros of Bs{u) deduced from analytic results in the large-/3o limit. Another example is the calculation 
of the renormalon residue at m = ^ in inclusive decays, which was first used as a constraint on Bs{u) in Ref. El]. 

For infrared and coUinear safe observables such as event-shape distributions, heavy-quark production cross section 
in e'^e~ annihilation, and heavy-meson decay spectra, one can solve the evolution equations 1)17(1 and ((25(1 and write 
a single factorization-scheme-invariant Sudakov factor. For example, in the case of e^e~ — > b -\- X the jet mass 
scale is Q^{1 — x) where Q is the e^e~ center-of-mass energy, while the soft scale associated with radiation off the b 
quark is m(l — x) where m is the b quark mass. The resulting Sudakov factor is 38]: 



Sud((5, m, N) — exp 



Cf 

l3o 



du 



T{u) 



A2 



Bsiu) 



T{~2u)T{N) 



-Bj{u) 



1 

2u 



2u) 

n-u)riN) 



(27) 



T{N-u) 

The corresponding anomalous dimensions Bs{u) and Bj{u) are now known to NNLO [53 . 15^ IssIIiOtIIiKtI i. A similar 
expression, with the same anomalous dimension functions^^ holds for inclusive B decay spectra, see Eq. (3.43) in 
Ref. Ei- 

It should be emphasized that Eq. ((27(1 is an exact all-order formula for the Sudakov factor. In particular, it is 
exactly renormalization-group invariant. In contrast with fixed-logarithmic-accuracy formulations, see e.g. there 
is no need to introduce any arbitrary renormalization or factorization scales. Naturally, uncertainties associated with 
unknown higher-order corrections to the anomalous dimensions Bg{u) and Bj{u) translate into uncertainties in the 
calculation of Sud(Q, m, N). By varying the functional form of Vj{u) and Vs{u) in Eqs. I(19() and ((26(1 . respectively, 
under the available constraints, one has a handle on the remaining uncertainty of the perturbative calculation. 

Although the DGE formulation and conventional fixed-logarithmic-accuracy approximations make use of the same 
perturbative calculation of the anomalous dimensions, the results are typically very different^^. These differences 



■^*The Sudakov factors of the heavy— quark distribution and fragmentation functions are the same to all orders |52| . 

■^^ Detailed numerical comparison between DGE with Principal Value prescription and fixed— logarithmic— accuracy approximations was 
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can be understood in detail by expanding the Borel function and examining the divergence of the emerging series |^ 
[2^. An important lesson from this comparison is that the order around which the renormalon divergence 
sets in varies significantly with N . Thus, there is no single optimal truncation order that could approximate the 
Borel sum. Moreover, when the asymptotic nature of the series sets in early, as occurs for example for inclusive 
decay spectra [s^. Issf. the logarithmic-accuracy criterion becomes completely useless: at that point the distinction 
between logarithmically-enhanced perturbative terms and power terms^^ is no more clear; separation must then be 
done without expansion. We note in passing that fixed-logarithmic-accuracy approximations have spurious Landau 
singularities in moment space which strongly restrict their range of applicability. Such singularities do not occur 
upon performing the integral in Eq. I|27|l . 

In order to uniquely define the perturbative Sudakov factor, and thus explicitly separate between perturbative and 
non-perturbative contributions at the power level, all renormalon singularities need to be regularized. The natural 
regularization is Principal Value integration. This regularization guaranties that the Sudakov factor, just like the 
physical moments, is a real- valued function of N , namely. 



Sud(P^)(Q,m,7V) = 



Sud(P^)(Q,m,7V*) 



(28) 



This implies that the corresponding pertubative spectrum, computed by Eq. ©, is real. 

Having defined the perturbative sum we are ready to consider power corrections. Power corrections on the soft 
scale Q/N appear in the exponentiation kernel H25|l and therefore generate a multiplicative exponential factor: 



where ^^^^^^ (Q) = cxp <( - >] ^(^) {N ~ 1){N ~ 2) . . . {N ~ k) ) , (29) 



with Zs 




Here non-perturbative corrections are introduced according to the ambiguity structure. Each renormalon pole in 
the evolution kernel (|25|l gives rise to a power correction with the corresponding A^-dependent residue. Zeros in the 
anomalous dimension Bs {u= |), as detailed in the third column in Tabled imply a missing renormalon ambiguity, 
so no corresponding power term is required. Note that the dimensionless parameters e^^ are expected to be of 0(1), 
so high powers in the sum arc inherently suppressed by 1/fc!, a direct consequence of the T{—2u) suppression factor 
in Eq. H25|l . This means that truncation of the sum in Eq. H29() is likely to be a good approximation up to high 
moments and a small number of non-perturbative parameters should suffice. 

It is important to emphasize that although S^'^^^^Q) sums up power corrections, its effect on the spectrum near 
threshold is not expected to be small. In general, the leading corrections become large for N ^ Q/A. Upon taking 
the inverse Mellin transformation the effect of S^^^^^\Q) can be recast as a convolution with a "shape function". 
Locally, this convolution generates an effect of 0(1); the most obvious example is a shift of the entire distribution 
by the k = 1 term in Eq. H29I) . Nevertheless, the formulation in Eq. H29|l is significantly more predictive than the 
conventional "shape function" approach, where the entire "soft" function Sn{Q,IJ-f) is being parametrized. The 
fundamental difference is that in Eq. (|29|) the first approximation to Sn{Q, fJ'p) is determined by a calculation and 
requires no non-perturbative input. Factorization scheme and scale dependence are completely avoided. Moreover, 
certain features (moments) of this distribution are protected by symmetries (Bg (u = §) =0) receiving no non- 
perturbative corrections. Finally, as explained above, the number of such non-perturbative corrections is limited by 
the factorial suppression inherent to high powers in Eq. H29(l . 



made in several examples, see e.g. Figs. 2 and 3 in Ref. |23for the case of the thrust distribution; Figs. 6 through 8 in Ref. ^3 ™ the 

context of heavy-quark fragmentation; and Fig 13 in Ref. I5i3l in the case of the triple differential rate in B > X^W decays. 

Consider for example contributions to the exponent in Eq. 1271 from u ^ ^. 
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Let us turn now to concrete examples. The fourth column in Table summarizes the power corrections appearing 
in the exponent in each class of inclusive distributions. In the following we discuss them one by one. We will focus on 
the comparison between the predicted power corrections and the ones observed in experiment wherever a dedicated 
DGE-based study was done. 

4.1. Drell-Yan Production 

As shown in Table ^ the case of Drell-Yan production near partonic threshold provides an interesting example 
where a u — ^ singularity, which could have appeared js^, 0| in l|^5|l is not realized once the detailed dynamics is 
taken into account :32j. Let us explain this conclusion using the DGE terminology established above [2^ : 

• The vanishing of Bg^{u — — 0; is understood to be a general property of this anomalous dimension, 
namely Bg"^{u = i) = in the full theory. This means that there is no m = ^ ambiguity in the corresponding 
evolution kernel. 

• The perturbative Sudakov exponent admits infrared safety at the 0{NA/Q) power level. Therefore, there is no 
reason to expect any 0{N K/Q) non-perturbative corrections. 

• The perturbative Sudakov exponent does have a -u = 1 ambiguity (see Tabled, indicating an O ((iVA/Q)^) 
power correction. More generally, TableQlshows that there are no odd power ambiguities, while all even powers 
are present. This renormalon structure leads to a specific parametrization of power corrections according to 
Eq. (|29|l . where the sum is over even values of k. In practice one should truncate this sum, restricting the 
number of parameters according to the kinematic range and the required accuracy. 

As already emphasized, the absence of a particular renormalon ambiguity does not imply the absence of corre- 
sponding non-peturbative effects, and therefore addressing the problem by alternative theoretical means is u seful . 
In the case of Drell-Yan, an alternative approach based on joint resummation [s^ leads to similar conclusions jll4| : 
there is no 0{NK/Q) correction, nor higher odd powers, while even power corrections are present. 

4.2. Event-shape distributions 

Event-shape distributions, being infrared and collinear safe jet observables 0|, provide a unique laboratory for 
resummation and power corrections. Consequently, these observables attracted much theoretical interest 
High-quality measurements of event-shape distributions in e^e~ annihilation have been performed over a wide 
kinematic range by several experiments, providing stringent tests of the theory. 

The basic motivation is to understand the effect of hadronization in a quantitative way. The starting point is 
a perturbative calculation of the distribution, where the variable is expressed in terms of on-shell quark and gluon 
momenta. The corresponding moments are finite owing to infrared and collinear safety. The computed spectrum is 
then compared with experimental measurements, where the variable is computed in terms of hadron momenta. 

One would like to understand the relation between the perturbative distribution and the measured one through 
power corrections. As usual, the main challenge is the description of the threshold region, the two-jet limit. This 
was the main motivation for much of the theoretical work on the subject and specifically, for developing the 

DGE approach [ii'23]. 

As shown in Tabled the single-dressed-gluon calculation predicts the absence of renormalons at all integer values 
of w; it therefore predicts only odd power corrections O {{N K/Q)^^ in the exponent. The ensuing phenomenology is 
the following: 

• The leading correction O (NA/Q), generates a global shift of the distribution and a, A/Q correction to average 
values of event-shape variables. This idea quickly and led to successful phenomenology: the shift 
can be described in terms of a single, approximately universal non-perturbative parameter in a large class of 
distributions '25'|. 
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Figure 7: Event-shape distributions: theory vs. data. The thrust (left) and the heavy-jet-mass (right) distributions, computed 
by DGE '2^ for center-of-mass energy of Q = Mz, are compared with LEP and SLD data. In the thrust case the three lines 
show: a Principal Value DGE calculation with no power corrections (dots); one which is shifted, where the shift parameter 
is determined is a global fit to data (dashed); and one that includes higher power corrections, determined in a fit (full line). 
These corrections are found to be consistent with the predication that the width does not change due to hadronization. In the 
heavy-jet-mass case no additional fits are performed: the dotted line is the result of a Principal Value DGE calculation and 
the full line includes hadronization corrections based on the fit to the thrust. 



• When approaching the threshold region, conventional Sudakov resummation accompanied by a shift of the 
distribution was proven insufficient. Here DGE [2^ [2^ made the difference. As shown in Fig. [7| the shape 
of the (Principal Value) DGE calculation, before any power corrections are included, is already very close to 
that of the data. Specifically, the width of the distribution is very similar, in accordance with the renormalon 
prediction that an O [{NA/Q}'^) correction would not occur. 



A dedicated DGE-based analysis in Ref. [23 has demonstrated that the measured hadronization effects are consis- 
tent with the renormalon predictions for power corrections. Successful fits to experimental data have been performed 
over the entire two-jet region and over a large range of energies. Moreover, as shown in Fig. \7\ the relation between 
hadronization corrections to different shape variables has been tested well into the threshold region. 

It is worthwhile noting that the value of af^(Mz) extracted from fits to the thrust distribution employing DGE [3| 
(or from the average thrust, employing renormalon resummation |20| ') is consistently lower than the world average 
value (by about 7%). The theoretical uncertainty (estimated [2^ as ±5%) is significantly larger than the experimental 
one (±1%), and it is dominated by uncertainty due to higher- orde r perturbative corrections. It is therefore important 
to revisit the comparison with data once 0(0;^) calculations |ll9| become available. 

An important comment is that in contrast with other distributions considered here, event-shape distributions are 
not completely inclusive Ii!tI|. For example, the thrust and the heavy-jet-mass distributions are defined based 
on the separation of space into two hemispheres. The single-dressed-gluon calculation of the Sudakov exponent, on 
which the power-correction analysis is based, assumes that the momentum of final-state particles associated with 
a single off-shell gluon is restricted to one hemisphere. This means in particular that certain correlations between 
the masses of the two hemispheres are neglected. The no-correlation assumption is also used when relating the 
thrust and the heavy-jet-mass distributions as shown in Fig.|7|above. The successful comparison with data indicates 
that these correlations are indeed small in the two-jet limit. 
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4.3. Heavy-quark fragmentation 



Consider next the fragmentation process in which an ofF-sheU heavy quark of mass m, produced at high energy 
0{Q), fragments into a heavy hadron. This process is infrared and coUinear safe. Therefore, the fragmentation 
function can be computed in perturbation theory, assuming that the final state is an on-shell heavy quark [stI Is^ 
I391 01 . This perturbative fragmentation function differs from the physical one by power corrections jsj E3 • 



XI 

\ 
b 
T3 



3 — 



~1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 T" 

solid: DGE with Principal Value 
dashes: NLL 
dots: no resummation 

Vq^ = 91.2 GeV 
m = 4.75 GeV 
A^^^ = 0.243 GeV 



0.0 




Figure 8: Heavy-quark fragmentation: perturbation theory vs. data. Comparison between perturbative calculations |39( of 
the normalized differential distribution in e^e~ — !• h -\- X with respect to xe = 2E/Q, and experimental data from Aleph. 
The dotted line is based on a simple NLO calculation with no soft-gluon resummation: only evolution logarithms ln{m^ /Q^) 
are resummed; the dashed line includes in addition soft-gluon resummation to NLL accuracy; the full line is the corresponding 
DGE result computed with the Principal Value prescription (the DGE Sudakov factor is quoted in Eq. 1271 above). 



Traditionally, heavy-quark fragmentation effects have been described in terms of "fragmentation functions" |63tl64 
a given functional form with one or more free parameters. Upon excluding the difficult threshold region where the 
fragmentation function peaks, such fits can indeed be performed. However, since there is no direct relation between 
these models and the QCD definition of the fragmentation function, the universality of the extracted parameters 
is doubtful. Within the threshold region such fragmentation models simply fail to bridge the gap between the 
resummed perturbative calculation and the data. This provided a strong incentive to apply the DGE approach to 
this problem [sglllot. 

As shown in Table ^ power-correction analysis of the Sudakov exponent in the large-/3o limit has definite predic- 
tions: 

• The leading infrared renormalon is located at u = i. Therefore, one expects a global shift of the distribution 
in Xe going from the partonic to the hadronic level. Comparing the data to the Principal value DGE result in 
Fig. |S1 one indeed observes such a shift. The shift parameter e^^^ in Eq. 129(1 can of course be determined by 
fitting the data, as shown in Fig. |51 below. 

• The u = 1 renormalon is absent. Therefore, one expects that the width of the spectrum would not be modified 
by non-perturbative fragmentation effects. Again, this is confirmed by the data js^; see Figs. |S| and El below. 

and above are present, and therefore higher-power corrections are 



• Higher renormalon ambiguities at u = 
expected. Their effect is restricted to high moments [3; 
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N Xj, 

Figure 9: Heavy-quark fragmentation: DGE-based fit to fieavy-meson production data from LEP. Left: results of a fit for 
A (i.e. Ois) and ef^ in Eq. (1291 (setting e^^ = for fc > 2) to Aleph moments A'^ = 2 to 6. Right: the corresponding curve in 
XE space, compared to the Aleph and Delphi data. 

Fits performed in Ref. [s^ confirm that Eq. H29|l above with a couple of leading non-perturbative power corrections, 
yields an excellent description of experimental data to very high moments. As shown in Fig.El a good description of 
data (and a consistent determination of ag) is achieved already with a single non-perturbative parameter, a plain shift 
of the PV-regularized DGE result. It should be emphasized that such power-correction phenomenology could not 
be developed starting with fixed-order calculations or fixed-logarithmic-accuracy Sudakov resummation: as shown 
in Fig. |S1 there is a qualitative difference between these results and the DGE one. 



4.4. Inclusive decay spectra 

As discussed in Sec. El the application of DGE to inclusive B decay spectra has been very successful. The most 
striking fact is that the on-shell approximation directly provides a good approximation of the spectrum. In other 
words, at present accuracy, power corrections on the soft scale O [{NA/Q)''^, the ones that distinguish between the 
quark distribution in a meson and that in an on-shell quark '5(J|, are negligible. 

To understand this issue, let us return once more to Table 

• As in the case of heavy-quark fragmentation (as well as in event-shape distributions) the leading infrared 
renormalon in the Sudakov exponent is located at u = ^. However, a unique property of inclusive decay spectra, 
is that this ambiguity cancels out exactly upon computing the resummed spectrum in physical kinematic 
variables [s^, EE 13 Ell • For example, for the radiative decay the resummed on-shell calculation yields: 



1 dT 



Ftnt dE~, 



2^ /-c+^oo ^ 



H{mb) JN{mb; ^f) x SNimb^fiF) + 0{1/N) 
^ ^ ^ 

Sud(mb,N) 



(30) 



Here S'Ar(mf,; /ii?) is the solution of Eq. I|25|) . the Sudakov factor describing the momentum distribution in an on- 
shell b quark. It has a leading ambiguity of the form exp{±(A^ — l)A/TOb}. The mass appearing i n the inve rse 
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UJ at 



Mellin transformation in H30() is the pole mass, which also has a leading renormalon ambiguity 
u = ^. Since the mass is raised to the power A'^ — 1, it generates an ambiguity of the form exp {t(A^ — l)A/mb} 
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in Eq. H30() . which cancels against the ambiguity of SNimi,; ^p). This cancellation was confirmed by explicit 
calculations of the corresponding renormalon residues in the large-/3o limit (Sec. 2.2.3 in Ref. [sojl'). 
Eq. I|3()(l refers specifically to the radiative decay, in which case the meson mass does not enter the calculation. 
A similar cancellation occurs in other inclusive distribution, see e.g. Eq. (4.4) in Ref. for the charmless 
semileptonic decay. There the ambiguity cancels upon expressing the partonic kinematic variables in terms of 
hadronic ones, a transformation that involves the mass difference between the meson and the quark. 

• As in the heavy-quark-fragmentation and event -shape-distribution examples considered above, the u = 1 
renormalon in Eq. (|25l) is absent. Therefore, one expects that the width of the spectrum would not be modified 
by non-perturbative Fermi-motion effects. Also here, comparison with data confirms the renormalon prediction: 
the photon-energy variance in Fig. computed as a function of the cut, is well described by the on-shell 
approximation. 

• Finally, there are renormalon ambiguities at u = | and above, which imply that non-perturbative power 
corrections due to Fermi motion should be included in Eq. H29I) with fc > 3. These should have an effect for 
high moments N ^ mt/A. Given the factorial suppression of high powers in Eq. (|29|l . a few non-perturbative 
parameters would be sufficient for any practical purpose, even if very precise data become available. 

The obvious conclusion is that having sufficient control on the perturbative calculation and the parametric form 
of power corrections, the paramctrization of the "shape function" and the associated uncertainty can be completely 
avoided. 

5. Conclusions 

The QCD calculation of inclusive cross sections and decay spectra near kinematic thresholds is challenging because 
of large perturbative and non-perturbative corrections. Conventional Sudakov resummation is limited to the range 
where non-perturbative corrections arc negligible. The DGE approach has proven effective in extending the range of 
applicability of perturbation theory well into the threshold region, where non-perturbative corrections are relevant. 
DGE makes maximal use of the known all-order structure of the Sudakov exponent and the inherent infrared safety 
of the observable. The main characteristics of this formalism are: 

• Factorization and exponentiation in moment space, going beyond the perturbative (logarithmic) level. 

• Resummation of running-coupling effects using Borcl summation, leading to rcnormalization-group invariance 
of the Sudakov exponent. 

• Power-like separation between the perturbative result of the exponent, computed as the Principal Value Borel 
sum and non-perturbative power corrections, whose parametric form is deduced from the renormalon residues. 

DGE has a high predictive power: 

• Important features of the distribution arc captured by resummed perturbation theory. 

• The perturbative result is modified by non-perturbative power corrections of a definite parametric form, see 
Eqs. (|20() and l|29() . The power-correction expansion in H29|l converges well even for high moments, and thus a 
small number of non-perturbative parameters is sufRcient. 

The main underlying assumption, whose validity cannot be addressed within perturbation theory, is that the dominant 
non-perturbative effects in the threshold region are indeed the ones identified by renormalon ambiguities. Here 
comparison with experimental data is essential. In all inclusive distributions considered so far, the observed pattern 
of power corrections is consistent with the renormalon prediction. Specifically, 
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• The potentially leading, u = i renormalon ambiguity in the "soft" Sudakov exponent Sn{Q', f^p) is of 
0(^{N — 1)A/Q). This corresponds to a global shift of the resummed distribution upon going from the partonic 
to the hadronic level. This shift has been firmly established and measured in event-shape distributions [j^. 
(see also 0,0,E3) in heavy-quark fragmentation [s^. In Drell-Yan or Higgs production the u = \ 
renormalon is absent due to a specific symmetry property of the corresponding anomalous dimension, which 
can be deduced [i^js^ from Eq. H26I) and Table. HI above (see an alternative approach in Ref. 34] ). In inclusive 
decay spectra the u = ^ ambiguity cancels out within the on-shell approximation, upon expressing the spectrum 
in terms of hadronic momenta (5(X i53t i55i | . 

• The potential u = 1 renormalon ambiguity does not occur in several cases, owing to symmetry properties of 
the Sudakov anomalous dimensions, which can be deduced from Eq. H26|l and Table. The corresponding 
non-perturbative correction would have induced an 0(1) change in the width of the spectra. Such effects 
have been excluded by comparison with data in event-shape distributions (see Ref. [23.E^ and Fig. [T] above), 
heavy-quark fragmentation (see Ref. [s^ and Fig. |H1 and 1^1 above) and inclusive decay spectra (see Ref. |54j 
and the variance plot in Fig. 13 above). 

The recent application of DGE to inclusive decay spectra has led to a fundamental change in the perception of 
the predictive power of resummed perturbation theory. For over a decade, it has been repeatedly stated that the 
calculation of the B — > Xs^ and B — > X^lv spectra is strictly beyond the limits of perturbation theory. The 
description of the spectra in the peak region has been based on "shape-function" phenomen olog y, \yhere an arbitrary 
functional form is introduced whose first few moments are constrained by data. Refs. have 
shown that resummed perturbation theory provides a good approximation to these spectra while genuinely non- 
perturbative effects (Eq. (|29|) with fc > 3) due to the "primordial" Fermi motion are not large. This advancement is 
already reflected in precision determination of mf^ and |Vub| using measurements of the B factories js^l- Moreover, 
in this framework there are high prospects for improving the determination of these parameters by incorporating 
higher-order (NNLO) calculations and by quantifying the relevant non-perturbative corrections. 
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